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SO ' ABSTRACT 

The scheme for accurate quantitative treatment of the radiation from a crystalline undulator in presence of 
the dechanncling and the photon attenuation is presented. The number of emitted photons and the brilliance 
'■^ ' of electromagnetic radiation generated by ultra-relativistic positrons channeling in a crystalline undulator are 
calculated for various crystals, positron energies and different bending parameters. It is demonstrated that with 
O ' the use of high-energy positron beams available at present in modern colliders it is possible to generate the 
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crystalline undulator radiation with energies from hundreds of keV up to tens of MeV region. The brilliance of 
the undulator radiation within this energy range is comparable to that of conventional light sources of the third 



c/2 ■ 

. generation but for much lower photon energies 
^ ■ Keywords: crystalline undulator, dechanneling, photon attenuation, brilliance 

ft ! 1. INTRODUCTION 

In this paper new results from the theory of electromagnetic radiation emitted by a bunch of ultra-relativistic 
positrons channeling through a periodically deformed crystal (a crystalline undulator) are reported. We formulate 
I , the approximation for effective analytical and numerical analysis of the characteristics of the undulator radiation 
with account for the influence of two main parasitic effects, the positron dechanncling and the photon attenuation. 
, The developed formalism is applied to calculate the number of the emitted photons and the brilliance of the 
radiation formed in crystalline undulators. 

' In a crystalline undulator there appears, in addition to a well-known channeling radiation, the radiation of 

an undulator type which is due to the periodic motion of channeling particles which follow the bending of the 
crystallographic planes. The parameters of the undulator radiation can be easily varied by changing the energy of 
beam particles and the parameters of crystal bending. The feasibility of this scheme was explicitly demonstrated 
for in Refs. 1, 2. In these papers as well as in the subsequent publications 3-8 the idea of this new type of radiation, 
the essential conditions and limitations which must be fulfilled to make possible the observation of the effect were 
formulated in an adequate form for the first time. A number of corresponding numerical results were presented 
to illustrate the developed theory. The importance of the ideas suggested and discussed in the cited papers has 



also been realized by other authors resulting in a significant increase of the number of publications in the field 
during last years 9-17 but, unfortunately, often without proper citation. 11-17 A detailed review of the results 
: obtained in this newly arisen field as well as a historical survey of the development of all principal ideas and 



related phenomena can be found in Ref. 18. 

The mechanism of the photon emission by means of a crystalline undulator is illustrated in Fig. 1. The 
(yz)-plane in the figure is a cross section of an initially linear crystal, and the z-axis represents the cross section 
of a midplane of two neighbouring non-deformed crystallographic planes (not drawn in the figure) spaced by the 
interplanar distance d. 

Further author information: A.V.S. is on leave from Ioffe Physical- Fechnical Institute, Russian Academy of Sciences, 
Polytechnicheskaya 26, St. Petersburg 194021, Russia. 
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Figure 1. Schematic representation of a crystalline undulator. Two black circles denote the nuclei belonging to two 
neighbouring crystallographic planes (separated by an interplanar distance d) which are periodically bent. The centerline 
of this channel (solid line) is described by a harmonic function y(z) — asm(2-Kz/\ u ). Its period A u and amplitude a satisfy 
the condition A u S> a. The dashed curve represents the trajectory of a projectile trapped in the channel. 

Under certain conditions the ultra-relativistic positrons will channel in the periodically bent channel. The 
trajectory of a particle contains two elements. Firstly there are channeling oscillations due to the action of the 
interplanar potential. Their typical frequency fl c h depends on the positron energy e and the parameters of the 
interplanar potential. Secondly there are oscillations related to the periodicity of the distorted midplane, - the 
undulator oscillations, whose frequency is uj = 2nc/X u . 

The spontaneous emission of photons is associated with both of these oscillations. The typical frequency 
of the channeling radiation is aj c h ~ 27 2 i7 c h where 7 = e/mc 2 is the relativistic Lorenz factor. The undulator 
oscillations give rise to the photons with frequency lo u w 4^ 2 luo/(2 + p 2 ) where p is the undulator parameter, 
p = 2nj(a/ X u ). 1,2 If strong inequality loq <C f2 c h is met than the frequencies of the channeling radiation and 
the undulator radiation are also well separated, oj <C w c h- In this case the characteristics of the undulator 
radiation are practically independent on the channeling oscillations but depend on the shape of the periodically 
bent midplane. 3,4 

There are essential features which distinguish a crystalline undulator from a conventional one based on the 
action of the periodic magnetic (or electric) field on the projectile. In the latter the beam of particles and 
the photon flux move in vacuum whereas in the proposed scheme they propagate through a crystalline media. 
Therefore, to prove that the crystalline undulator is feasible, it is necessary to analyze the influence of the 
interaction of both beams with the crystal constituents. On the basis of such analysis one can formulate the 
conditions which must be met and define the ranges of parameters (which include e, a, A u and also the crystal 
length L u and the photon energy hw) within which all the criteria are fulfilled. In full this analysis was carried 
out very recently and the feasibility of the crystalline undulator was demonstrated in an adequate form for the 
first time in Refs. 1, 2 and in Refs. 3-8. 

For further referencing let us briefly mention the conditions which must be met in a crystalline undulator. 

A stable planar channeling of an ultra-relativistic positron in a periodically bent crystal occurs if the maximum 
centrifugal force, F c {, is less than the maximal force due to the interplanar field, F lnt . Notating the ratio F c {/Fi nt 
as C one formulates this condition as follows 1 ' 2 ' 19 : 



C = (2tt) 2 
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< 1. 



There are two essentially different regimes of the radiation formation in a periodically bent crystals. They are 
defined by the magnitude of the ratio a/d. In the case of low amplitudes, a/d <C 1, the characteristic frequencies 
of the channeling radiation and the undulator radiation become compatible uj u ~ ui c h- This results in the loss of 



the monochromaticity of the radiation, since the channeling radiation is essentially non-monochromatic due to 
noticeable deviations of the interplanar potential from a harmonic form. Additionally, in this case the intensity 
of undulator radiation is small compared with that of the channeling radiation. 3 ' 4 

On the contrary, in the limit a 3> d not only the characteristic frequencies are well separated, w u /o> c h ~ 
C d/a <C 1, but also the undulator radiation intensity is higher than the intensity of the channeling radiation. 1 ' 2 ' 4 
As a result, if one is only interested in the spectral distribution of the undulator radiation, one may disregard the 
channeling oscillations and assume that the projectile moves along the centerline of the bent channel. Therefore, 
the criterion which is imposed on the relative magnitudes of d, a and A u is as follows 

d<a<A u . (2) 

The second inequality ensures that the crystal is deformed elastically, and its structure and symmetry arc not 
affected by the deformation. 

The term 'undulator' implies that the number of undulator periods, N u , is large. Only in this limit does 
the radiation formed during the passage of a bunch of relativistic particles through a periodic system bear the 
features of undulator radiation (narrow, well-separated peaks in spectral-angular distribution) rather than those 
of synchrotron radiation. Hence, the following strong inequality, which entangles the period A u and the length 
of a crystal L u must be met in the crystalline undulator 1,2 : 

iV u = ^ » 1. (3) 



The coherence of the radiation, emitted in the crystalline undulator, takes place if the energy of the channeling 
particle does not change noticeably during with the penetration distance. For ultra-relativistic projectiles the 
main source of energy losses are the radiative losses. Therefore, it is important to establish the range of energies 
for which the parameters of undulator radiation formed in a perfect periodic crystalline structure are stable. In 
Ref. 3 a comprehensive quantitative analysis of the radiative loss of energy, Ae, due to the channeling and the 
undulator radiation was carried out. It was established that the relative radiative losses Ae/e become large if 
the initial energy of the positron bunch is e > 10 Gev. For lower energies of positrons 

e < 10 GeV, (4) 

the radiative losses are small, Ae < O.Ole. 

As was pointed out Rcfs. 1, 2, 5-7 two phenomena, the dechanncling effect and the photon attenuation, lead 
to severe limitation on the length of a crystalline undulator. 

If the dechanneling effect is neglected, one may unrestrictedly increase the intensity of the undulator radiation 
by considering larger 7V u -values. In reality, random scattering of the channeling particle by the electrons and 
nuclei of the crystal leads to a gradual increase of the particle energy associated with the transverse oscillations 
in the channel. As a result, the transverse energy at some distance from the entrance point exceeds the depth 
of the interplanar potential well, and the particle leaves the channel. The mean penetration distance covered 
by a channeling particle is called the dechanneling length. For given crystal and projectile the dechanncling 
length Ld = Ld(e, C) depends on the energy e and on the parameter C (see (1)). To calculate the dechanneling 
length one can either apply the diffusion theory to describe the multiple scattering 20,21 or carry out a computer 
simulation of the scattering process of the projectile from the crystal constituents. 5 ' 22 Alternatively, to estimate 
id(e, C) one can use the approximate formulae. 2, 22 For an ultra-relativistic positron the dechanncling length in 
straight channels (i.e. C = 0) in various crystals lies within the interval Ld(e, 0) (cm) w (0.05 . . . 0.08) e (GeV), 18 
i.e. does not exceed several millimeters at GeV energies of a positron. For a periodically bent channel the 
dechanncling length decreases as C grows following, approximately, the law L&{e, C) « (1 — C) 2 id(e, 0). 5,22 

The propagation of photons emitted in a crystalline undulator is strongly influenced by the atomic and the 
nuclear photoeffects, the coherent and incoherent scattering on electrons and nuclei, the electron-positron pair 
production. All these processes lead to the decrease in the intensity of the photon flux as it propagates through 
the crystal. A quantitative parameter, which accounts for all these effects and defines the scale within which the 



intensity of a photon flux decreases by a factor of e, is called the attenuation length, L a = L & {lo). It is related 
to the mass attenuation coefficient /j,(u>) as L a (u>) = l/fi{u). 23,24 The coefficients n(u>) are tabulated for all 
elements and for a wide range of photon frequencies. 23 The magnitude of L a (ui) depends on u and on the type 
of the constituent atoms. For low-Z crystals (e.g., diamond) the magnitude of L a {ui) exceeds that for a tungsten 
crystal taken for the same ui by a factor of 10 1 . . . 10 3 . In the case of a diamond crystal the value of L a (ui) varies 
from 1CP 2 cm at Tilo w 5 keV up to several cm for hu > 10 5 eV. 

The simplest way to account for the dechanneling and the attenuation is to consider the case when the crystal 
length satisfies the condition L u < min [L&(e, C), L a (u>)], and to assume that within the chosen L u scale neither 
the number of channeled particles nor the flux of emitted photons do decrease. Such approach was utilized in 
most of the papers devoted to the crystalline undulator problem. More consistent treatment of the dechanneling 
process and of its influence on the parameters of the undulator radiation was carried out in Rcfs. 5, 6, where a 
simple analytic expression for the spectral- angular distribution was derived which contains, as a parameter, the 
dechanneling length. 

In the present work we make another step in developing the theory of the crystalline undulator. The following 
problems are solved and discussed below in the paper. 

(a) We propose the scheme for accurate quantitative treatment of the radiation from a crystalline undulator 
in presence of the dechanneling and the photon attenuation (Sect. 3). As a result, we evaluate an analytic 
expression for the spectral- angular distribution and the number of emitted photons which contains, as 
parameters, three quantities L u , L&(e,C), and L a (uj). 

(b) We demonstrate that for given type of the crystal and crystallographic plane, and for given values of e, 
a, A u and w there exists an optimal length of the crystal which ensures the largest number of the emitted 
photons. 

(c) Using (a) and (b) we carry out the calculation of the number of emitted photons and the brilliance of the 
crystalline undulator radiation (Sect. 4). The calculations, which account for the conditions (l)-(4), are 
performed for several crystals and by using the parameters of positron bunches used in modern colliders. 24 

Prior to the discussion of the radiation formed in a crystalline undulator, in Sect. 2 we briefly summarize the 
results from the general theory of undulator radiation (see, e.g. Refs. 25-27). We use the term 'ideal undulator' 
to indicate that the propagation of positrons and photons occurs in vacuum. 



2. CHARACTERISTICS OF RADIATION FORMED IN AN IDEAL UNDULATOR 

The spectral-angular distribution of the energy E emitted by an ultra-relativistic projectile in a planar undulator 
can be written in the following form 



d 3 E 
Tidw dfl 



S(u,6,<p) D N ^(fj). 



(5) 



Here 9 <C 1 and <j> are the emission angles with respect to the undulator axis, dSl = QdOdip is the solid angle of 
the emission. The function S(u>,0,ip), which does not depend on the undulator length, is given by 



S{u,6,<p) = £t -4^1 {p 2 \h\ 2 + 7 2 e 2 \I \ 2 - 2 P1 9 cos tp Re (I*h)} , 
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(6) 
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Here a w 1/137, uio = 2irc/\ u , p is the undulator parameter and the parameter r\ is given by 



The factor Dn u (t)) on the right-hand side of (5) is defined as follows 



DnM 



sin7V u 7r77 x 2 
sin nfj 



(9) 



where fj = rj — n and n is a positive integer such that n — 1/2 < r\ < n + 1/2. 

For N u 3> 1 the function Z)jv u (??) has a sharp and powerful maximum in the point fj = 0, where -Djv u (0) = ./V 2 . 
The width of the peak A?7 U is equal to 1/N U . This behaviour of Djv u (??) results in a peculiar form of the spectral- 
angular distribution of undulator radiation which clearly distinguishes it from other types of electromagnetic 
radiation formed by a charge moving in external fields. Namely, for each value of the emission angle 9 the 
spectral distribution consists of a set of narrow and equally spaced peaks (harmonics). The peak intensity is 
proportional to TV 2 . This factor reflects the constructive interference of radiation emitted from each of the 
undulator periods and is typical for any system which contains N u coherent emitters. 

The values w„ of the harmonics frequencies follow from the condition that parameter r\ becomes an integer 
(this corresponds to fj = 0). In particular, in the case of the forward emission, 8 = 0, the harmonics frequencies 
are defined from the relation 



1 u n I l + P_ 

27 2 wo V 2 



(10) 



For 9 — the integrals (7) can be evaluated analytically and the spectral-angular distribution calculated for 
uo = uj n (for n = 1,3,5...) acquires the form 25, 28 : 



d 3 E 
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(l+p 2 /2) 2 



J n-1 (z) — J n + 1 (z) 



(11) 



where z = np 2 / (4 + 2p 2 ) and J v {z) is the Bessel function. 

The finite width of the central peak of Dat u (?7) defines the emission cone Af2„ and the bandwidth Aut n /uj n 
of the nth harmonic. Using Afj u = l/N u and accounting for (10) one derives 



7T l+ P 2 /2 

7 2 nN u 



(12) 



Formulae (11)-(12) allow one to calculate the number of photons AiV Wn of energy u> 
Ato n /2 emitted by a beam particle within the cone AQ n : 

d 3 E 



AN^ 



hdudQ 



a r> AuJn at n r \ AcJr > 



u>„ - Aw tl /2,w n + 



(13) 



where Q n (p) = 4z [J(„-i)/ 2 (z) - J^ n+1 )/ 2 (z)] . 

Let us introduce two other quantities which characterize the radiation formed in an undulator and are closely 
related to the number of the emitted photons, but also take into account the properties of the beam of ultra- 
relativistic particles. These quantities are the flux and the brilliance (see, e.g. Ref. 27). 

The flux F n describes the number of photons per second of the nth harmonic emitted in the cone AQ n 
and in a given bandwidth. A quantitative definition of this quantity, measured in (photons/s/0.1%Bw) (the 
abbreviation 'BW stands for the bandwidth Au) n /w n ), is given by the following formula 28 : 



F n = 



T7S37T^7 \ ~ = W 3 ™N u Q n (p)- = 1.431 x 10 14 N u Q n (p) I [A] , 
10 3 (Aw„/w„) e e 



(14) 



where / is the electric current of the beam. In the latter expression I is measured in Amperes. 



The general definition of brilliance of the photon source of a finite size is given in terms of the number of 
photons of energy huj emitted in the cone Ail per unit time interval, unit source area, unit solid angle and per 
bandwidth. 27 To calculate this quantity is it necessary to know the beam sizes a x , a y and angular divergencies 
<fr x , cj>y in two perpendicular directions, as well as the divergency angle of the radiation and the 'size' of the 
photon beam. The brilliance of undulator radiation can be related to the flux F n as follows 28 : 

B " = <sfc/ (15) 



Here e XtV = + (J xy y^ 2 + <fr xy are the total omittance of the photon source in the x and y directions, with 

(fr n = ^/Af2„/27r being the angular width of the nth harmonic and a n — \ n /Air(fr n is the 'apparent' source size 
taken in the diffraction limit. 29 

To obtain brilliance in the units (photons/s/mrad 2 /mm 2 /0.1%Bw) the quantities a x . v and a n must be 
measured in millimeters whereas the angular variables (fr x . y and <fr n - in milliradians. 



3. CHARACTERISTICS OF RADIATION FORMED IN A CRYSTALLINE 

UNDULATOR 

In an ideal undulator the beam of particles and the emitted photons propagate in vacuum. In a crystalline 
undulator, due to the interactions with crystal atoms, the particles can dechannel, and thus be lost for further 
motion through the undulator. Additionally, the photons emitted inside the crystal can be absorbed or scattered 
while making their way out from the crystal. Therefore, it is necessary to account for the processes of dechanneling 
and photon attenuation. In what follows we carry out the qualitative analysis of the influence of these two 
processes on the characteristics of the radiation formed in a crystalline undulator. 

3.1. Spectral-angular distribution in presence of the dechanneling and attenuation 

Let the crystal length, the amplitude and period of bending and the energy e satisfy the conditions (l)-(4). 

A positron, which enters the crystal at small incident angle with respect to the curved crystallographic plane, 
penetrates through the crystal following the bending of its channel. However, due to random scattering by the 
electrons and nuclei of the crystal the energy of the transverse oscillations of the positron in the channel increases, 
and finally the particle leaves the channel, becoming lost for the crystalline undulator. Although the rigorous 
treatment of the dechanneling process cannot be implemented by analytical means only, it is possible to develop 
a model approach based on the assumption that the probability w(z) for a particle to penetrate at a distance z 
along the undulator axis (z <E [0, L u ]) can be described by the exponential decay law 

w(z) = cxp(-z/i d ) . (16) 

In intermediate formulae when referring to the dechanneling length we omit its arguments e and C. 

With the effect of dechanneling taken into account the spectral-angular distribution of the radiated energy 
per one particle can be written as follows: 

d 3 E _ d 3 £(-> | d 3 £(+> 
Ti&uj dfl hduj dD, hduj dfl 

The first term is the contribution to from all the processes in which the particle dechannels somewhere inside 
the crystal. To calculate this term one notices that the quantity Ld^dzexp (—z/L^) defines the probability of 
a particle to channel through the distance z and then dechannel within the interval dz. Such a particle emits 
the radiation which corresponds to the undulator of the length z and the number of periods z/X u . Therefore 

d 3 £(-> _ f L dz_ _ z/Ld d?E^ att \z) 
hdujdfl J Ld hdudfl 



(18) 
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Figure 2. Illustration of the photon attenuation in a crystalline undulator. A photon (the long-dashed line), emitted 
within the jth period of the undulator of the length z, can be absorbed (or scattered) in the part of crystal of thickness 
L u — jA u on its way to a distant detection point Ro (Ro 3> L u ). 



where d 3 'E^ att ^ (z) / 'fidw dfi denotes the spectral-angular distribution from the undulator of the length z. The 
superscript '(att)' indicates that to calculate this quantity one has to account for the photon attenuation. 

The second term on the right-hand side of (17) is due the process when the projectile channels through the 
whole length L u . Its probability is given by the factor exp (—L u /Ld). Therefore one can write 

^ (+) = e - W , d d3 ^ (aU) (^) (m 



If the photon attenuation is neglected, then to calculate (19) one uses (5) instead of d 3 E^ att \L u )/hdujdn. 
The integral in (18) is also evaluated with the help of (5) where one substitutes iV u with z/X u . Such approach 
was applied Ref. 5 with the only difference that in the cited paper to calculate (18) and (19) we used the 
discrete probabilities instead of the continuous distribution function (16). The use of the latter implies that the 
dcchanncling effect is small over the scale of one undulator period and, therefore, » A u . 

Now let us turn to the derivation of the quantity d 3 i^ att ) (z) /Tidoj dil which is the spectral-angular distribution 
of radiation formed in the undulator of the length z < L u in presence of the attenuation. In the intermediate 
formulae we assume that the ratio N z = z/ A u is an integer number which corresponds to the number of periods 
in this undulator. In the final formula this limitation will be omitted. Throughout the text the notations L u 
and N u are reserved for the length of the crystal and the number of undulator periods within L u . 

As mentioned above, if one neglects the photon attenuation effect, the distribution d 3 E^ att \z) / hduj dQ. is 
described by (5) where one substitutes N u with N z . The only quantity in (5) which depends on the number 
of undulator periods is the factor D^ z (fi) defined in (9). This factor appears in the formula for spectral- 
angular distribution as a result of squaring the modulus of a coherent sum of the amplitudes of electromagnetic 
waves emitted from spatially different but similar parts of the undulator. In more detail, Dm z (v) IS given by 

Dnz(v) — Sjii cxp^ifci? ~ 2i7rryj^ . The argument (kR Q — 2TTrjj) (k = lj/c is the wavenumber) stands for the 
phase of the electromagnetic wave emitted within the jth period of the undulator and detected at some distant 
point Rq from the undulator. It is assumed that the quantities L u , z and Rq satisfy the relations: z < L u <C Rq. 



In a crystalline undulator a photon emitted within the jth period in the direction of the point R can be 
absorbed within the distance L u — j\ u while propagating through the crystal, see Fig. 2. To account for this 
possibility one can assume that the wavenumber becomes complex, k — > uj/c + i/i/2. The quantity /i = /i(a->) 
defines the attenuation length L a (uj) = ^^(u) within which the photon flux is reduced by a factor of e. For a 
complex k the phase factor e lkR ° , which in an ideal undulator is the same for all periods j = 1 . . . N z , is replaced 



with e lkR °e ^ Lu i Xu \ and a proper expression for Dn u is 
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1 + c mA u - 2c^ A "/ 2 cos(27T?7) ' 
The spectral-angular distribution of radiation in presence of the photon attenuation acquires the form 



(20) 



hdco dCl 



S(L0,8, i p)D ( £ t) (rj) 



(21) 



In the limit \i — > (i.e., when there is no attenuation) the factor D^ tt \fj) becomes equal to Dpj z (fj) from (9), 
and the right-hand side of (21) reduces to that of Eq. (5). 

To derive the explicit expression for the spectral-angular distribution of the radiated energy from a crystalline 
undulator one uses (21) in (17)-(19). Let us note here, that although the expression (20) was obtained for the 
case when the ratio z/X u is an integer, its use in the integral from (18) can be justified by the above mentioned 
conditions, that the undulator period A u is small compared to the dechanneling length L d , the attenuation length, 
L a , and the length of crystal, L u . Hence, the relative error, which appears when one uses (21) in (18), is small, 
being of the order of magnitude A u / min{Ld, L a , L u } <C 1. Carrying out the integration one represents the total 
spectral-angular distribution (17) of radiation formed in the crystalline undulator in the form similar to (5) 



d 3 L 
hdu) dSl 



S{w,0,<p)VnM: 



(22) 



where the function S(uj 7 8 7 ip), defined by (6), does not depend on L u , Ld and L a . These parameters enter the 
factor T>N u {fj) which is given by the expression: 
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—2 | cos 



-K d _ 
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K d (2«d - Ka) 2 + 40 2 
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d (2« d -«a) 2 +40 2 



where the following notations are used: 



K d 



Lu 
L' 



27rry7V u . 



(23) 



(24) 



Despite a cumbersome form of the right-hand side of (23) its main features can be easily understood. Firstly, 
we notice that if the dechanneling is neglected, id — ► oo (or Kd — > 0), the function T>N u (fi) reproduces D^ tt \fj) 
from (20). In another limit Kd = K a = (i.e., no attenuation and dechanneling) eq. (23) reduces to the definition 
of the factor D/v u (?)) which characterizes the ideal undulator. In the case when only the attenuation effect is 
neglected the limit of X>jv u (r?) can also be easily evaluated. In either of these cases the main maximum of 2?at u (?)) 
is located in the point fj = 0, i.e. when the parameter rj reduces to an integer, and, therefore, the harmonics 
frequencies are still defined by (10). The maximum value T>n u {0) can be presented as follows: 



IW a (0)=4JV d 2 



2 1 



+ 



(l-x)(2-x) x(l-x) x(2-x) 



where the quantity iVd = Ld/A u stands for the number of undulator periods within Ld, and the ratio 



x = 



Kd 



Ld 
L a 



(25) 



(26) 



does not depend on the crystal length L u . 



The width of the central peak A77, which in the case of an ideal undulator equals to 1/N U , is increased 
due to the photon attenuation and the dcchanncling. Formally, the additional widths are due to the factors 
V( K a + sin 2 71-77) an d l/((2«d — K a) 2 + 4</> 2 ) which enter (23). The widths associated with these factors are, 
respectively, A771 = K a /(27V u 7r) and Af]2 = |2«d — K a |/(2AT u 7r). Thus, the total width of the peak is: 



Ar? = ^u- 2 + (A771) 2 + (Ar) 2 ) 2 = j-Jl + ( " a ffi + (27) 

The additional widths lead to the enlargement of the solid angle AQ n of the emission cone in the forward 
direction. In accordance with (27) one derives 



AQ - = 7^rV 1+ ^^-n^- < 28 > 

The formulae for the number of photons AAf Un emitted in the cone Att n as well the corresponding flux of 
radiation T n one derives similarly to how it was done in Sect. 2 for an ideal undulator. The result is: 

AA/L„ - naN eS (x, K d ) Q n (p) ^ (29) 

T n = 1.431 x 10 14 N cS (x, «d) Q„(p) / [A] . (30) 

The difference between these equations and formulae (13) and (14) is that the number of undulator periods N u , 
met in the latter, is substituted with the effective number of periods, N c ff(x, Kd), which is defined as follows: 



v nM .A , .2 (^-i) 2 + i 



N eS (x,K d ) = ^^^l + Kp ^ = N d f(x,K d ) (31) 



f(x,K d ) = — 
K d 



+ 
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(1 - x) (2 - a;) x(l - x) x(2 - x) 
We use these equations in the subsequent section to define the optimal length of a crystalline undulator. 
3.2. Optimal length of a crystalline undulator 

In the case of an ideal undulator one can, in principle, increase infinitely the length of the undulator. This will 
result in the increase of the number of photons, the photon flux, and the brilliance since they are proportional 
to the number of periods. The limitations on the values of L u and N u are mainly of a technological nature. 

The situation is different for a crystalline undulator, where the number of channeling particles and the number 
of photons which can emerge from the crystal decrease with the growth of L u . It is seen from (32) that if L u — > 00 
then the parameters Kd = L u /Ld and XKd — L u /L a also become infinitely large, and the effective number of 
periods goes to zero leading to AM 0Jn ,J r n — > 0. This is quite natural result, since in the limit L u 3> Ld all 
particles leave the channeling mode and, thus, do not undulate in the most part of the crystal, whereas all 
emitted photons are absorbed inside the crystal if L u 3> L & . Another formal (and physically trivial) fact, which 
follows from (31) and (32), is that N e s(x, Kd) = also for a zero-length undulator, when L u = 0. Vanishing of a 
positively defined quantity N e ff(x, Kd) at two extreme boundaries suggests that there exists the length L(x) for 
which the effective number of periods (taken for fixed values of L a , Ld and A u ) attains the maximum. 

To define the value of L{x) or, what is equivalent, the quantity Kd{x) = L(x)/Ld, one carries out the derivative 
of f(x, Kd) with respect to Kd and equalizes it to zero. The analysis of the resulting equation shows that for 
each value of x = Ld/L & > there is only one root R d . Hence, the equation defines, in an inexplicit form, a 
single-valued function Rd(x) = L(x)/Ld which ensures the maximum of N e ff(x, Kd) for given L a , Ld and A u . 

It is important to note that the crystal length enters Eqs. (29)-(30) only via the ratio Kd- All other quantities, 
met in these formulae as well as in (31) and (32), are independent on the length of the crystal. Therefore, the 
quantity L(x) ensures the highest values of AA/L n and T n for the radiation formed in the crystalline undulator. 
In this sense L(x) can be called the optimal length of the undulator which corresponds to a given value of x. 
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Figure 3. Dependences Rd(x) = L(x)/Ld and f x {x,Kd(x)) = N e ff(x,Kd(x))/Nd on x = Ld/L a . 

The dependences of Rd(x) — L(x)/Ld and of the ratio f(x, Rd(x)) = N c g(x, Rd(x))/Nd on x are presented in 
Fig. 3. For a given crystalline structure, the dechanneling length Ld is uniquely defined by the energy e and the 
parameters of bending a and A u . On the other hand, the attenuation length L a is the function of w. Therefore, 
fixing e, a, A u and u> one calculates x — L&/ L a and, then, using the dashed curve in the figure finds the optimal 
length of the crystalline undulator L(x) which accounts for the dechanneling effect and the photon attenuation. 
Simultaneously, from the solid curve one finds the effective number of the undulator periods N D ff(x, Rd{x)) which 
defines the number of emitted photons, the flux and the brilliance of radiation. 

4. NUMERICAL RESULTS 

From (29) follows, that to find the number of photons AA/L>„ emitted in a crystalline one has to calculate two 
factors. The factor Q n (p) (see (13)) depends on the harmonic number n and on the undulator parameter p 
which, in turn, is defined by the values of e, a and A u through the relation p = 2n^a/X u . The second factor, 
N e g(x, Kd) depends on L u , A u , = Ld(s,C) and L a = L a (w). It was explained in Sect. 3.2 that once the 
quantities A u , Ld(e,C) and L a (uo) are known the length of the crystal can be fixed by the condition L u = L 
which results in the maximum values of N e g(x, Kd) and T n with respect to L u . The numerical data presented 
below in this section was obtained for the optimal length of undulator. 

Therefore, to calculate AJ\f Un one fixes, in addition to the crystallographic plane, the values of n, e, a and 
A u (the three latter are subject to the conditions (l)-(4)) which uniquely define the quantities p, C, w n , L a (ui). 
However, there is some uncertainty with respect to the magnitude of the dechanneling length. This uncertainty is 
not intrinsic to the case of a periodically bent crystal but rather reflects the stochastic nature of the interaction 
of a channeling particle with crystal constituents. As mentioned above to calculate the dechanneling length 
one can apply the diffusion theory to describe the multiple scattering or carry out numerical simulations of the 
scattering process. Alternatively, one can use model-dependent analytic expressions for Ld(s, C). In the present 
paper we utilize the approach, presented in Ref. 2, and approximate Ld{e, C) with 

L d (s,C) = (l-C) 2 Ld(e,0), L A (e, 0) = ^ £ . (33) 




Here where ro = 2.8 x 10 - cm is the electron classical radius, rac 2 = 0.511 MeV is the electron rest energy, (Ztf 
is the Thomas-Fermi radius of the crystal atom. The parameter C is defined by Eq. (1). The quantity Ld(e,0) 
stands for the dcchanncling length of a positron in a straight crystal. 5 ' 22 The quantity A = In v /2^mc 2 /I—23 /24, 
with / denoting the (average) ionization potential of the crystal atom, is the Coulomb logarithm characterizing 
the ionization losses of an ultra-rclativistic particle in amorphous media. For a quick estimation of Ld{s,0) 
(in cm) one can re- write the right-hand side of the second equation from (33) as 2aTF^£/A, with axF and d 
measured in A and e in GeV. The values of qtf and d , are presented in Tabic 1. 



Table 1. Parameters d, o,tf and I7m ax for different crystals and channels. 





C (111) 


Si (111) 


Ge (111) 


W (110) 


d(A) 


1.54 


2.35 


2.45 


2.24 


a T F (A) 


0.258 


0.194 


0.148 


0.112 


EC X (GeV/cm) 


9.23 


8.58 


17.5 


57.4 



To calculate the brilliance of a crystalline undulator (which one obtains by using (30) in (15)) it is necessary to 
specify the parameters of a positron bunch, which are the current J, the beam sizes a XtV and angular divergencies 
<px,y We used the parameters of the positron beams from several modern high-energy e~e + colliders. These 
parameters are summarized in Table 2. The data on e, cr XtV , I, Af and I (which is an average beam current) are 
taken from Ref. 24. The beam divergencies <j> XtV were calculated using the data on the transverse emittance (not 
presented in the table) and the beam size <J x ,y The peak current /, which is defined as the electric current of a 
single bunch, was calculated as / (A) w 48A/"/ 1 with / in cm. 

Table 2. Positron energy e, bunch length I, number of particles per bunch Af, beam sizes ff X;y , beam divergencies 4> x , y , 
and a positron peak current I for several modern high-energy e~e + colliders. 24 





DA$NE 


VEPP-2000 


BEPC-II 


PEP-II 


KEKB 


CERS-C 




(Frascati) 


(Russia) 


(China) 


(SLAC) 


(KEK) 


(Cornell) 


s (GcV) 


0.700 


1.0 


1.9-2.1 


2.5-4 


3.5 


6 


I (cm) 


1-2 


4 


1.3 


1 


0.65 


1.2 


Af (units 10 iu ) 


3-9 


16 


4.8 


6.7 


7.3 


1.15 


a x (mm) 


0.800 


0.125 


0.380 


0.157 


0.110 


0.300 


<j y (mm) 


0.0048 


0.125 


0.0057 


0.0047 


0.0024 


0.0057 


(j) x (mrad) 


0.375 


2 


0.379 


0.153 


0.164 


0.500 


4> y (mrad) 


0.208 


2 


0.544 


0.319 


0.417 


0.439 


1(A) 


144-216 


192 


177 


322 


539 


46 



The results of our calculations arc presented in Figs. 4 and 5. The choice of the crystals was motivated by 
the fact that C, Si, Ge and W crystals are frequently used in channeling experiments (see, e.g., Ref. 22). An 
additional reason is that for a given photon frequency the magnitude of L a (uj) rapidly decreases with the growth 
of atomic number of the constituent atoms. Therefore, by comparing the results obtained for different crystals 
one can investigate the influence of the photon attenuation on the formation of the radiation in a crystalline 
undulator. 

Graphs 1 (a) ' in Fig. 4 present the dependence of the maximal number of emitted photons of the first harmonic 
(n = 1) per bandwidth Acji/wi and per a positron versus the ratio a/d. The curves were calculated for the 
positron energies indicated in Table 2 (for BEPC-II and PEP-II colliders we used the values e = 2 GeV and e = 3 
GeV, respectively). For each crystal and for each e value the dependences ( AJV" Wl /BW) max were obtained as 
follows. There are two independent variables, A u and a, which, (for fixed crystal, energy and harmonic number n) 
define all other quantities on the right-hand side of (29). For practical purposes it is more convenient to chose the 



ratio a/d > 1 and the parameter C < 1 (see Eq. (1)) as the independent variables. Then, for each pair (a/d, C) one 
finds A u , p = 2n'ya/X u , Qi(p), the dechanneling length L&{e,C) and the number of periods = Ld(s,C)/X u , 
the fundamental harmonic frequency u)\ (see Eq. (10)) and the attenuation length L a (wi), and the value of 
N c ff(x,R,d(x)) which corresponds to the optimal undulator length calculated for x = La(e,C)/L a (oji) (see (31)- 
(32) and Fig. 3 and Sect. 3.2). As a result, one finds the magnitude of AA/^/BW. Finally, scanning through all 
(a/d,C) values one determines the highest possible value of the number of photons per BW, (AjV u)1 /BW) max , 
as a function of a/d. Having done this one also finds the dependence uii — ui(a/d) (graphs '(b)' in Fig. 4) as 
well all other characteristics of the undulator as functions of a/d. 

Let us briefly discuss the behaviour of obtained dependences. Firstly, as it is seen from the graphs (a), for 
a fixed amplitude a the quantity (AJ\f Ul /BW) max is an increasing function of a positron energy e. This feature 
becomes clear if one analyzes the e dependence of the product Q\{p) f{x, Rd{x)) which defines the number of 
emitted photons (see Eqs. (29) and (31)). All three factors are increasing functions of energy (although it is not 
too obvious for f(x,Rd(x))). 

Another feature of the curves (AA/" u , 1 /BW) max is that they are decreasing function of a/d in the region 
a/d > 1. To a great extent this is a consequence of the photon attenuation in the crystal. Indeed, as the ratio 
a/d increases the undulator period A u increases too, in order to maintain the inequality C« 1 (see Eq. (1)). 
Larger values of A u results in lowering of the emitted photon energy (see Eq. (10) and the graphs (b) in Fig. 4) 
and, consequently, to the decrease of the attenuation length, L a {uj). This, in turn, leads to the increase of the 
ratio x — Ld/L a which defines the magnitude of f(x,Rd(x)). This factor, as it is seen from Fig. 3, rapidly falls 
off for x > 0.1, and this feature manifests itself in the dependence of (AA/" Wl /BW) max on a/d. In the case of 
crystals consisting of heavy atoms the dependence acquires additional features, which are due to the fact that 
the ionization potentials, Io, of the inner atomic subshclls of such atoms lie within the energy range 1 . . . 100 
keV. The photons with the energy just above the threshold are absorbed much more efficiently than those with 
the lower energies. As a result, the dependence of L a (uj) in the vicinity of the threshold becomes a saw-like. For 
u> < I the attenuation length noticeably (up to the order of magnitude) exceeds L a (oj) for w > I . This effect 
results in the irregularities of the dependence (AAd/BW) max on a/d, which in Figs. 3 are mostly pronounced 
for diamond and tungsten crystals. 

In the opposite limit, when a/d <C 1 the number of the emitted photons goes to zero. This tendency, which 
is seen explicitly for all the curves (but the CERS-C one) in the case of C and Si crystals, is also clear and is due 
to the fact that the case a — corresponds to the linear crystal, i.e. the absence of the crystalline undulator. 

In Fig. 5 we present the peak brilliance of the crystalline undulators based on different crystals (as indicated) 
and calculated using the parameters of the positron beams from Table 2. The data refer to the emission in the 
first and the third harmonics in the forward direction. It is seen that in contrast to the number of the emitted 
photons which is the same, by the order of magnitude, for all colliders, the magnitudes of the peak brilliance for 
different beams differ by orders of magnitude. To the largest extent this is due to the quality of the beam, which 
includes, apart from the beam current /, its size and angular divergency, see. Eqs. (14) and (15). For all crystals 
and over the whole range of photon energies the product e x e y of the photon source emittances is the smallest for 
the KEKB collider (labeled as '5' in the graphs in Fig. 5). As a result, this beam, which does not lead to the 
highest values of (AA/" wl /BW) max , ensures the largest peak brilliance of the crystalline undulator radiation. The 

peak brilliance for the KEKB positron beam is on the level of (4 . . . 20) x 10 22 ^photons/s/mrad 2 /mm 2 /0.1%Bw) 

for the photon energies within 1 ... 10 MeV range. These values can be compared with the peak brilliance of 
the light sources of the third generation. 30 The peak brilliance on the level 10 21 . . . 10 23 in the 100 keV range 
of photon energies by means of the undulators based on the action of magnetic field is planned to be achieved 
within several projects. 31-33 The data from Fig. 5) demonstrate that it is feasible to produce the radiation of 
the same level of brilliance but for much higher energies by means of crystalline undulators. 

5. CONCLUSION 



Theoretical investigations show that it is entirely realistic to use a crystalline undulator for generating sponta- 
neous radiation in a wide range of photon energies. The parameters of such an undulator, being subject to the 
restrictions mentioned in Sect. 1, can be easily tuned by varying the parameters of the bending, the positron 



energy and by choosing different channels. The large range of energies available in modern colliders together 
with the wide range preparation of periodically bent crystalline structures allow one to generate the crystalline 
undulator radiation with energies from hundreds of keV up to tens of MeV region. The brilliance of the undulator 
radiation within this energy range is comparable to that of conventional light sources of the third generation but 
for much lower photon energies. 

The experimental efforts are needed for the verification of numerous theoretical predictions. Such efforts will 
certainly make this field of endeavor even more fascinating than as it is already and will possibly lead to the 
practical development of a new type of tunable and monochromatic radiation sources. 
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Figure 4. Graphs (a): the maximal number of photons of the first harmonic (n = 1) per a bandwidth Acji/oji and 
per a positron as a function of the ratio a/d calculated for the positron energies in various colliders (see Table 2) as 
indicated. Graphs (b): the corresponding values of the fundamental harmonic energy (see Eq. (10) with n = 1). See 
also explanations in the text. Each vertical pair of the graphs (a) and (b) correspond to the positron channeling in the 
particular periodically bent channel as indicated in the graphs (a). The legend refers to all graphs in the figure. 
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Figure 5. Peak brilliance of the undulator radiation in the forward direction calculated for four channels as indicated in 
each graph. The solid curves correspond to the radiation in the fundamental harmonic n = 1, the dashed curves refer to 
n — 3. In each graph the enumerated sets of the solid and the dashed curves correspond to the parameters of the positron 
beams in different colliders (see Table 2). 1: DA$NE, 2: VEPP-2000, 3: BEPC-II, 4: PEP-II, 5: KEKB, 6: CERS-C. 



